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ABSTRACT 



The supernova remnant (SNR) RX J0852-4622 (Vela Jr., G266.6-1.2) is one of the most important SNRs for investigating the acceler- 
ation of multi-TeV particles and the origin of Galactic cosmic rays because of its strong synchrotron X-ray and TeV y-ray emission, 
which show a shell-like morphology similar to each other. Using the XMM-Newton archival data consisting of multiple pointing ob- 
servations of the northwestern rim of the remnant, we investigate the spatial properties of the nonthermal X-ray emission as a function 
of distance from an outer shock wave. All X-ray spectra are well reproduced by an absorbed power-law model above 2 keV. It is found 
that the spectra show gradual softening from a photon index T = 2.56 in the rim region to F = 2.96 in the interior region. We show that 
this radial profile can be interpreted as a gradual decrease of the cutoff energy of the electron spectrum due to synchrotron cooling. 
By using a simple spectral evolution model that includes continuous synchrotron losses, the spectral softening can be reproduced with 
the magnetic field strength in the post-shock flow to less than several tens of fiG. If this is a typical magnetic field in the SNR shell, 
y-ray emission would be accounted for by inverse Compton scattering of high-energy electrons that also produce the synchrotron 
X-ray emission. Future hard X-ray imaging observations with Nustar and ASTRO-H and TeV y-ray observations with the Cherenkov 
Telescope Array (CTA) will allow to us to explore other possible explanations of the systematic softening of the X-ray spectra. 
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Supernova remnants (SNRs) have been considered as major 
cosmic -ray (CR) accelerators below the knee energy of ~ 
10'^ eV. The measured energy density of CRs can be explained 
if ~ 10% of each supernova kinetic energy is transferred to the 
accelerated particles. A plausible mechanism for this particle ac- 
celeration is diffusive shock acceleratio n (DSA; e.g.. lBelllll978l; 
lprurvll983l:rBlandford &. Eichler<1987l) . In the DSA theory, par- 
ticles are accelerated at the shock, resulting in a power-law dis- 
tribution, which can account for the Galactic CR spectrum using 
a reasonable diffusion coefficient in the Galaxy. However, many 
unsolved problems are still remaining, such as the acceleration 
efficiency, the maximum energy of the accelerated particles, and 
so on. Observations of synchrotron X-ray emission have been 
playing important roles in stimulating the development of the 
DSA theory. For example, magnetic field amplification at the 
shocks of yo ung SNRs ha s been inferred from X-ray observa- 
tions (e.g.. lUchivama et al.i 12007) and it is now considered as 
an int egral part of the DSA theory (e.g., lEllison & Vladimirovl 
l2008h . 

RX J0852.0-4622 (also called G226.2-1.2 or Vela Jr) is a 
young SNR with ~ 2° of an angular size in the line of sight to 
the Vela SNR. RX J0852. 0-4622 was originally discover ed in the 
data of the ROSAT All-Sky Survey ("Aschenbach" 1998V Based 
on the ROSA T Asia, whose bandp ass was limited to soft X-rays 
below 2 keV, lAschenbachl (119981) argued that the emission from 
the SNR can be explained either by a hot thermal model with a 



temperature of ~2.5 keV or by a power-law model with a photon 
index of F ~ 2.6. ASCA observations, with an imaging capability 
of up to ~ 1 keV, showed that the X-ray emission from the bright 
rims of RX J0852. 0-4622 is dominated by a nonthermal emis- 
sion characterized by a power-law with a photon index F ~ 2.6 
(iTsunemi et al.ll2000l: [Slane et al.ll200ii) . The nonthermal X-ray 
emission is presumably of synchrotron origin. The Chandra im- 
age of the northwestern (NW) rim revealed shar p filamentary 
structures similar to those discovered in SN1006 (iBamba et al.l 
2005). 

The CANGAROO collaboration claimed the detection of 
TeV y-rays from the direction that coincides with the peak of 
the X-ray emission in the NW rim (Katagiri et_aL,2QQ5). Then, 
the H.E.S.S. telescopes, with their highly sensitive stereoscopic 
observations, have detected spatially extended TeV y-ray emis- 
sion, which sho ws a good spati al correlation with the nonthermal 
X-ray shells ( Aharonian et al. ^'2007). The H.E.S.S. spectrum of 
RX J0852.0-4622 can be well fitted with a power-law of F ~ 2. 1 . 
Recently, the Large Area Telescope (LAT) on board the Fermi 
Gamma- ray Space Telescope has detected y-ray emission from 
RX J0852.0-4622 in an energy band of 1-300 GeV (Tanaka et al] 
1201 II) . The LAT spectrum is described by a power-law with a 
photon index F ^ 1.85 + 0.06, which is substantially harder than 
the H.E.S.S. spectrum. The origin of the y-ray emission has been 
actively debated; the GeV-TeV y-rays can be explained either 
by inverse Compton (IC) scattering of high-energy electrons, or 
by ;7ro-decay y-rays produced in interact ions of accelerated pro - 
tons (and nuclei) with interstellar gas (lAharonian et al ] l2007h . 
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Fig. 1. EPIC-MOS image of the northwestern rim of RX J0852.0 -4622 in an energy range of 0.5-10 IceV (left) and the scheme to search for the 
radial profile of the synchrotron X-ray emission (right). The angular size of the summed image is larger than XMM-Newtons FOV. Numbers in 
figure indicate region IDs. A contamination source in region 8 was removed from the spectral analysis. 



SNR RX J0852. 0-4622 is one of the most interesting objects for 
studying CR acceleration at SNR shocks. 

The age and di stance of RX J0852. 0-4622 are still debated. 
Ilvudin etal.l(ll998h estimated an age of ~ 680 yr and the dis- 
tance of ~ 200 pc based on the flux of **Ti lines detected with 
COMPTEL on board the Compton Gamma-Ray Observatory. 
iTsunemi et al.l (12000 ) estimated a similar age of between 630 
and 970 yr based on observations of Ca X-ray lines with ASCA. 
These estimates imply that RX J0852. 0-462 2 belongs to the su- 
pernovae that are closest to Earth. However, ISlane et aP (1200 l b 
argued that the column density for the X-ray spectrum of the 
SNR is higher than that for the Vela SNR and that the distance to 
RX J( )852.0-4622 should be m uch larger, 1-2 kpc. On the other 
hand, iMoriguchi et al ] (1200 Ih observed the molecular distribu- 
tion using '^CO (J=l-0) emission measured with the millimeter 
and submillimeter telescope NANTEN. They estimated the up- 
per limit on the distance to be ~ 1 kpc. Recently, a new estimate 
of the age and distance has been reporte d from an expansion 
measurement of the NW rim of the SNR (iKatsuda et al. 2008); 
1700^300 yr and 0.75 kpc for the age and distance. Assuming 
a high shock speed of 3000 km-s these values seems probable 
for the SNR. 



2. Data and analysis 

2.1. Data reduction 

We used XMM-Newton archival data of 14 pointing observations 
with slightly different aiming points. We analyzed data from 
the European Photon Imaging Camera (EPIC), which consists of 
two MOS and one PN CCD arrays. All observations were per- 
formed in the PrimeFuUWindow mode, either with the medium 
or the thin filter The difference between the medium and the 
thin filter is excluded from the following analysis since we only 
used a high-energy band of 2-10 keV for spectral fitting. Reduc- 
tion and analysis of the data were performed following the stan- 
dard procedures using the SAS 11.0.0 software package. Since 
the fluxes measured by the EPIC PN instrument show systemat- 
ically lower values (by 10-15%) in comparison with the fluxes 
measured by the EPIC MOS, we used only the MOS data. The 



total exposure time amounts to 436 ks for EPIC-MOS 1. The 
summary of the archive data is shown in Table 1 . 

2.2. Spectral analysis 

To investigate the spatial properties of the nonthermal emission 
in the NW rim, we first combined different pointing images from 
EPIC-MOS 1 and MOS2 with the emosaic tool. Although an ex- 
pa nsion of the shock fron t over a time span of 6.5 yr was reported 
bv'K atsudaeTai] ( |2008|) . we did not consider this because the 
possible displacement (~ 0.1') due to the remnant expansion is 
quite small. Figure 1 shows the combined image in an energy 
range of 0.5-10 keV without exposure correction. Sharp edges 
that represent outer shock fronts can be clearly seen in the im- 
age. There is also a fainter shock structure located inside the 
outer boundary (see region ID 6 in Fig. 1). The X-ray emission 
arises from the interior past of the remnant, even ~ 10' from the 
rim. 

The source photons were accumulated from separate rect- 
angular regions with an angular width of 0'.67 in the region 1- 
4 and r.34 in the region 5-9 as shown in Fig. 1 (right). We 
co-added source photons of both detectors (MOS 1 and 2) to in- 
crease the statistics. Since thermal emission from the Vela SNR 
is dom inant below 2 keV , we only used 2-10 keV for spectral 
fitting (Hiraga et al. 2005). All spectra were well reproduced by 
an absorbed power-law model. Figure 2 shows the background- 
subtracted spectrum of each rectangular region. The absorption 
column density was fixed at A^h = 0.58 x 10^^ cm"-, which 
was determined from region 1 . Figure 3 (top and middle pan- 
els) shows the fitting results of the photon index and the flux 
integrated from 2 keV to 10 keV. The photon index shows grad- 
ual softening from 2.56 on the rim to 2.96 in the interior region. 
Though radial profiles of the surface brightness have been stud- 
ied in various shell-type SNRs, a gradual softening of the photon 
index like this has not been seen before. It provides a new clue 
as to how high energy electrons evolve in a shock downstream. 
The results are summarized in Table 2 (left column). 

The nonthermal X-ray spectrum in RX J0852. 0-4622 can be 
interpreted as synchrotron emission produced by high-energy 
electrons at a cutoff region of their energy spectiTim. The ob- 
served systematic change of the photon index in the radial direc- 
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Fig. 2. Spectrum of each rectangular region with the best-fit absorbed power-law model. Parameters are given in the left column of Table 2. 



tion (Table 2, left column and Fig. 3, top panel) is most likely 
caused by the change of the cutoff position in the synchrotron 
spectrum. Therefore, to investigate this idea, we fit the same 
data by another model. We employed the following function to 
fit the X-ray spectruniJ 



F{e)^€-^ ■exp[-i-f\ 
So 



(1) 



where F is the photon index and eo is the cutoff energy. In this 
approach, our aim is to determine eo from the observed varia- 
tion of F. Therefore the cutoff energy eo was treated as a free 
parameter, while the index F was fixed. This functional form 
can be regarded roughly as the synchrotron radiation produced 
by electrons with an energy distribution of 



N{E) oc E-P ■ exp [-(- — )], 



(2) 



where fimax is the maximum energy of electrons. Applying a 
(J-functional approximation to the synchrotron spectrum of a 
single-energy electron, we obtain Eq. (1) from Eq. (2) with a 
relation of F = (p + l)/2. Moreover, the cuto ff energy of the 
synchrotron spectrum eo is related to /Smax as (iAharonianll2000l: 
Ulchiyama et al. 2003) 



eo - 5.3 ■ ( 



B 



:)■( 



lOfiG' MOOTeV 



(3) 



In some young SNRs, a shock-acceleration spectrum with 
p ^ 2.2 has been inferred from recent multiwavelength data. 
Hence, one should expect F ^ 1.6, provided that the power-law 
part of the electron spectrum does not suffer from deformation 
due to energy-dependent losses. On the other hand, if the cool- 
ing time of the electrons responsible for the synchrotron X-ray 
emission is shorter than the source age, the volume-integrated 



' The high-energy cutoff of shock-accelerated electrons is due to the 
balance between acceleration and synchrotron radiat ion losses at the 
shock front if the magnetic field is strong. Zirakas hvili & AharonianI 
(|2007) have derived the energy spectrum of synchrotron radiation under 
such circumstances. To allow for other scenarios for the formation of 
the spectral cutoff, we employed Eq. (1). 



spectrum is expected to be characterized by a steeper power- law 
part of F ~ 2. We adopted two cases of F = 1.6 and 2.0. The 
best-fit cutoff energies are presented in Table 2 (middle and right 
columns) and Figure 3 (bottom panel). Quoted errors are at the 
Icr confidence level. We obtained good fits (x^, ji) for each re- 
gion. 
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Fig. 3. Variation of the photon index (top panel), the flux integrated 
from 2 keV to 10 keV (middle panel), and the cutofi" energy eo obtained 
from an exponential cutofi' power-law model with F = 1.6 (bottom 
panel). 
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Table 1. Summary of the archive data. 



Obs. ID 


Instrument Mode 


Filter 


Obs. Date 


MOSl GTI (ks) 


112870301 


PrimeFullWindow 


Medium 


2001 Apr 25 


30.4 


137550901 


PrimeFullWindow 


Medium 


2001 Dec. 10 


34.9 


153750701 


PrimeFullWindow 


Medium 


2002 May 21 


22.3 


156960101 


PrimeFullWindow 


Medium 


2002 Nov. 06 


13.6 


159760101 


PrimeFullWindow 


Medium 


2003 Jun. 22 


17.5 


159760201 


PrimeFullWindow 


Medium 


2004 Nov. 01 


22.4 


159760301 


PrimeFullWindow 


Thin 


2005 Nov. 01 


33.4 


162360101 


PrimeFullWindow 


Medium 


2003 Nov. 10 


23.9 


162360601 


PrimeFullWindow 


Medium 


2003 Nov. 10 


2.1 


162363101 


PrimeFullWindow 


Medium 


2003 Dec. 24 


18.7 


412990101 


PrimeFullWindow 


Thin 


2006 Nov. 09 


58.9 


412990201 


PrimeFullWindow 


Thin 


2007 Oct. 24 


61.2 


412990501 


PrimeFullWindow 


Thin 


2008 Nov. 12 


55.9 


412990601 


PrimeFullWindow 


Thin 


2009 Oct. 25 


40.7 



Table 2. Summary of the spectral fittng. 



power-law 



exp. power-law (F = 2.0) 



60 [keV] xfjyf 



exp. power-law (F = 1.6) 



60 [keV] 



Reg. ID A^H [-lO'Vcm 



Flux (2-10 keV) 



0.58 ±0.09 



2.56 + 0.02 
2.63 ± 0.02 
2.68 + 0.02 
2.74 ± 0.03 
2.84 ± 0.02 
2.91 ±0.03 
2.99 ± 0.05 

2.95 ± 0.09 

2.96 ± 0.22 



2.74 


±0.01 


1.08 (176) 


3.28 


±0.01 


1.12 (197) 


2.66 


±0.01 


1.03 (167) 


2.10 


±0.01 


0.93 (119) 


1.18 


± 0.01 


1.19(59) 


0.80 


± 0.01 


0.89 (59) 


0.47 


± 0.01 


1.14(58) 


0.29 


± 0.01 


1.2(30) 


0.16 


± 0.01 


0.90 (32) 



0.42 ± 0.08 3.59 ± 0.27 1.03 (176) 

2.84 ±0.18 1.13(197) 

2.43 ±0.17 1.06(167) 

2.10 ±0.20 0.98 (119) 

1.80 ±0.18 1.21(59) 

1.37 ±0.14 0.82(59) 

1.04 ±0.19 1.20(58) 

1.00 ±0.41 1.24(30) 

0.96± 0.63 0.97 (32) 



0.14 ±0.03 



1.36d 

1.14: 
1.01 : 
0.92: 
0.86: 
0.72 : 
0.66: 
0.55 : 
0.56: 



0.06 1.06 (176) 
0.03 1.09 (197) 
:0.03 1.02 (167) 
0.06 0.92 (119) 



0.06 
0.06 
0.07 
0.10 
0.12 



1.12(59) 
0.80 (59) 
1.17 (58) 
1.22 (30) 
0.98 (32) 




Angular distance [aronnin] 

Fig. 4. Spatial dependence of the cutoff energy 6o obtained from the 
exponential cutoff power-law model of F = 1.6 with the time-dependent 
spectral evolution model. 



3. Discussion 

3.1. Interpretation of variation in spectral shape 

The fitting results do not show any significant difference except 
for the absorption column density between simple power-law 
and exponential cutoff models because of a limited energy band. 
However, it is still hard to interpret that the electrons are formed 
in different power-law shapes with region-dependent accelera- 
tion efficiencies behind the shock front. This assumption comes 
from the widely accepted particle-acceleration theory, i.e., the 



diffusive shock acceleration or the first-order Fermi accelera- 
tion, which specifies the effective acceleration site at the shock 
front with a power-law electron spectrum with an index of ~ 2. 
Therefore the different photon indices can be more naturally in- 
terpreted as the presence of spectral cutoffs associated with the 
unavoidable cutoff in the acceleration spectrum of parent elec- 
trons. A solid evidence of such cutoff energy has been reported 
in a young SNR RX J1713.7-3946 (on the order of 1000 yr) with 
the Suzaku broadband spectr um of more than two decades in the 
energy band, i.e., ~ 40 keV (iTanaka et al.ll2008h . Although the 
FoV of the Suzaku H XD is not sm all enough to allow detailed 
imaging spectroscopv. lTanaka et a l. (2008)'s study suggests that 
a spectral cutoff is a common feature for the synchrotron spec- 
trum, which shows a typical photon index of 2-3 below 10 keV of 
shell-type SNRs. We note that the obtained absorption column 
density also indirectly suggests t he presence of the e nergy cutoff 
in the spectrum of the Vela Jr. ISIane et aP (1200 ih pointed out 
that the column density for the power-law is significantly higher 
than that f or Vela SNR, e.g., ~ 10^° cm'^ dBocchino etalJI 19991: 
lAschenba ch et al.lll995h . If we assume the distance to the Vela 
Jr to be 0.75 kpc, the lower column density with the exponential 
cutoff model, i.e., 1.4 x 10^' cm"^, seems physically preferable 
in comparison with a distance to the Vela SNR of ~0.25 kpc 
( ICha et ai]|1999 ). We thus interpret the different photon indices 
as energy cutoffs of the parent electron spectra in the following 
discussion. 
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Fig. 5. Radial profiles of the observed X-ray flux (top panel) and 
approximation models (bottom panel), proposed by Petruk et al. (201 1). 
The horizontal axis indicates the normalized distance from the outer 
shock wave. 



3.2. Synchrotron cooling effect on the cutoff energy 

The variation of the cutoff energy reflects the time evolution of 
the electron spectrum. Once a particle distribution has been in- 
jected at the shock, it evolves downstream, suffering from ra- 
diative losses and adiabatic expansion losses. We first ignored 
the adiabatic loss and focused on the synchrotron cooling effect 
on the cutoff energies to estimate an upper limit of the magnetic 
field strength. We consider the following equation to describe 
a time evolution of electrons in a certain volume element only 
subject to synchrotron losses: 



= ^[b(E)N(E)] + Q(E,t), 
at oE 



(4) 



where f denotes time elapsed since the shock, b{E) - aE^ is 
the energy loss rate of the electrons, namely, b{E) - -(^), and 
Q{E) is a source term that describes the rate of injection of elec- 
trons and their injection spectrum into the source region. If we 
assume the injection of electrons with a power-law energy spec- 
trum at f = without subsequent injection of electrons, we can 
write the electron spectrum as Q{E) - KE^f6{t), where 6{t) is 
the Dirac delta function. It is straightforwar d to show that the 
solution of Eq. (4) is given as (lLongairlll994h 



NiE, t) = QqE-p ■ (1 - atE)''-^ ■ 0(- - E), 

at 



(5) 



evolution of the electron energy by integrating the energy loss 
function, i.e., / = / a ■ dt, and then. 



imax = (at + Ej) \ 



(6) 



where Eo = E^^ySt - 0). The radial evolution of eo can be cal- 
culated from £'max(0 with a transformation of f into x, which is 
the distance from the shock to the volume element. We assume 
that the shock velocity has a time dependence of V oc r^^/^ j-^j^g 
Sedov-Taylor evolution), where T denotes the age of the rem- 
nant. The coefficient of the synchrotron energy loss rate a is 
given as 



b{E) = aE^ 

= 6.6x lOVfi^ [eV/s] 
B 2 

" ^-^^ ^TOT^G^ '^ToOTeV 



^ f [keV/s], 



(7) 



where y is the Lorentz factor of the electron, and B is magnetic 
field strength. Defining k - a/B^ , Eq. (6) can be written as 



E^,Ax)^[kB^-t + E^Y\ 



(8) 



where B is the magnetic field strength (assumed to be constant in 
t, i.e., B is uniform downstream, as well as T , i.e., B is uniform 
upstream of the forward shock). From Eq. (3) we can obtain the 
cutoff photon energy as 



£o(x)ocB-[kB^-t + Ef^^]~^. 



(9) 



We adopt a cuiTent shock velocity of V = 30 00 km-s ' 
and distance of D = 0.75 kpc ('Katsu da etani2008h . The co- 
effcient £0 is chosen so that the spatial curves pass the initial 
data point. Then, the magnetic field strength controls the spa- 
tial variation of the cutoff energy. The time-dependence of the 
shock velocity, i.e., the Sedov-Taylor evolution, is included in 
the model. We used Vi-ei = (l/4)y as a relative speed between 
the shock front and post-shock gas to transform t into x. Figure 4 
shows the spatial dependence of the cutoff energy with the time- 
dependent spectral evolution model. The magnetic field strength 
of B ~ lOyuG can reproduce the decrease of eo as a function of 
X under the assumption that the field strength does not change 
both in time and space. We note that a projection effect, which is 
neglected in our analysis, dilutes the spatial variation of eo^ and 
consequently could increase the required magnetic field strength 
somewhat. 



3.3. Adiabatic and inverse Compton scattering effects on the 
cutoff energy 

The electrons would significantly suffer from the adiabatic 
losses. In this case, the required magnetic field strength be- 
comes lower than the above estimated value, and then the inverse 
Compton scattering effect also needs to be considered. The adi- 
abatic loss rate for the electrons due to expansion of the volume 
is given as (Longair, 1994) 



where is a step function. The factor (1 - atEy^^ represents a 
cutoff shape at E^ut = (of) ' due to synchrotron losses. We may 
interpret iSmax of Eq. (2) as E^m. However, since the injection 
spectrum also has a cutoflF in reality, this interpretation would 
not be always appropriate. 

If we assume that the downstream evolution of /imax is sim- 
ply determined by synchrotron losses, we can obtain the time 



(-) 

V ''adiabatic 



1 .da 
~~R'^'dt' 



E, 



(10) 



where a is the fluid coordinate along the radius of the remnant 
R. Here we again assume that an electron population has been 
confined within a region of certain volume and injected at the 
shock in a fluid form. The energetic electron fluid does work 
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adiabatically as it expands, and consequently loses its internal 
energy. 

As for the inverse Compton scattering, the seed photon fields 
include the cosmic microwave background (CMB) and Galactic 
infrared and optical fields in the interstellar medium. In gen- 
eral, the photon fields produced locally by the SNR itself would 
not contribute to the target soft photons and Compton scattering 
of CMB photons should dominate the others. In the Thomson 
regime, the energy-loss rate via inverse Compton scattering of 
the CMB field is written as 



(^)lC = -gCt^TySVf/cMB, (11) 

where cr-j- is the Thomson cross-section, p is the velocity of an 
energetic electron in units of the speed of light c, and Ucmb - 
0.26 eV cm"^ is the energy density of CMB photons. We initially 
assumed an input electron spectrum as an exponential power-law 
shape with p - 2.2 and numerically calculated the electron spec- 
trum after losing energy due to the synchrotron, adiabatic, and 
inverse Compton scattering effects. The cutoff energies were ob- 
tained by fitting the synchrotron spectra produced by the contin- 
uously energy-modulated electron spectra. We assumed a uni- 
form interstellar magnetic field strength of ~ 1 /iG. As is shown 
in Fig. 4, the adiabatic effect dominates over the synchrotron 
cooling and the cutoff variation can be explained without radi- 
ation losses. In the adiabatic loss-dominant case, a site of the 
magnetic field amplification might be focused on the bright fil- 
ament region and the field strength of the inner regions is sim- 
ilar to a normal interstellar magnetic field. The magnetic field 
strength of the filament is estimated as ~ 6 jiG fr om the X- 
ray/T eV observations and the leptonic model ( Aha ronian et al.l 
l2007h . This implies that the magnetic field amplification near 
the shock occurs only moderately in the SNR. However, Eq. (10) 
might be a too simplified assumption because the electron popu- 
lation might not be suitable as an ideal expanding fluid for over 
1 kyr, e.g., at a post shock flow speed on the order of 3000/4 
km-s"', electrons fall behind the shock 1' or 6.74 x 10'^ km in 
~ 300 yr. Thus we conclude that the magnetic field strength near 
the filament is estimated to be between 6-10 /iG. 

3.4. Application of MHD-based model to tine radial profile 

To confirm the consistency of the above estimated magnetic field 
strength with an MHD- based calculation, we applied the analyt- 
ical model proposed by Petruk et aP (1201 ih to the obtained radial 
profile. This model accounts for the projection and other MHD 
effects, such as a variation of the magnetic field strength in the 
post-shock region, and calculates synchrotron images of SNRs 
of the Sedov-Taylor phase in a uniform interstellar medium. Fig- 
ure 5 shows the comparison between the observed radial profile 
of the X-ray flux and the model curves. We assumed the cut- 
off energy near the shock front to be eo = 4 keV and the shock 
obliquity angle ©o = 0°, which is defined as the angle between 
the magnetic field and the normal to the shock. For other pa- 
rameters, we used typical values as follows: k-^^ - \, - 8.25, 
o- = 4, r = 5/3, /7 = 0, ^ = 0, r = 1000 yr, ekev = 2 keV 
(see Petruk et al. 2009, 2011 for the definitions of the model 
parameters). 

We calculated radial profiles for three sets of ( B, fimax )'■ (i) 
B = 5 yuG and E^^„ = 122.8 TeV, (ii) B = 10 /zG and E^^^ = 
87 TeV, and (iii) B = 20 juG and ^.^x = 61.4 TeV. Each set 
gives eo = 4 keV using Eq. (3). The angular radius of the SNR 
was adopted to be /? = 1°. The radial profile of the X-ray flux 



is broadly consistent with the models with B - 5-20 //G, which 
roughly agrees with the magnetic field used above to reproduce 
the spatial variation of eg. 

If the y-ray emission detected by Fermi and H.E.S.S. is pre- 
dominantly caused by IC scattering, the magnetic field in the 
y-ray production region should be ~ 10/iG because of the mea- 
sured ratio of the TeV y-ray and X-ray energy fluxes, u)y (1-10 
TeV)/wx (2-10 KeV) ~ 2 (Aharonian et al. 2005). The magnetic 
field required to explain the radial profile would favor the lep- 
tonic model. However, this is not conclusive because the obser- 
vational constraints are not sufficient to explore other possibili- 
ties that may account for the variation of ey. More detailed mul- 
tiwavelength observations are necessary to place more definitive 
constraints on the magnetic field and to elucidate the y-ray emis- 
sion mechanism. 
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